Inelastic neutron scattering study of phonon density of states in nanostructured 

Sii_^Ge^. thermoelectrics 



Chetan Dhital/ D. L. Abernathy,^ Gaohua Zhu/ Zhifeng Ren,^ D. Broido,^ and Stephen D. Wilson^-j^] 

^Department of Physics, Boston College, Chestnut Hill, Massachusetts 024 67, USA 
^Quantum Condensed Matter Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6393, USA 

Inelastic neutron scattering measurements are utilized to explore relative changes in the gener- 
alized phonon density of states of nanocrystalline Sii_a:Gea; thermoelectric materials prepared via 
ball milling and hot pressing techniques. Dynamic signatures of Ge-clustering can be inferred from 
the data by referencing the resulting spectra to a density functional theoretical model assuming ho- 
mogeneous alloying via the virtual crystal approximation. Comparisons are also presented between 
as-milled Si nanopowder and bulk, polycrystalline Si where a preferential low energy enhancement 
and lifetime broadening of the phonon density of states appear in the nanopowder. Negligible 
differences are however observed between the phonon spectra of bulk Si and hot pressed, nanostruc- 
tured Si samples suggesting that changes to the single phonon dynamics above 4 meV play only a 
secondary role in the modified heat conduction of this compound. 

PACS numbers: 63.22.-m, 63.20.-e, 63.22.Kn, 65.80.-g 



I. INTRODUCTION 

Recent advances in materials synthesis techniques have 
enabled grain boundary and strain engineering in a num- 
ber of bulk, functional, nanocomposite alloys resulting 
in dramatic changes in their physical propertieJ^^. One 
prominent example of this has been the ability to engi- 
neer nanostructured thermoelectric composites with sub- 
stantially reduced thermal conductivities yet which re- 
tain reasonably robust conduction pathways, resulting 
in an en hancement in the thermoelectric figure of merit 
ZlH2llHini, While the enhanced efficiency of these ther- 
moelectric nanocomposites continues to be optimized,^ 
remarkably little is understood concerning the mecha- 
nism governing their thermal conductivity — the primary 
factor in their enhanced performance relative to their tra- 
ditional bulk counterparts. 

One of the simplest and most striking examples of 
an enhanced ZT in a bulk theromoelectric nanocompos- 
ite was recently reported in a series of Sii-a^Gca; alloys 
prepared via high energy ball milling and hot pressing 
technique^. A substantial drop in the thermal conduc- 
tivity of P-doped, nanostructured. Si can be achieved 
relative to bulk Si, and this gain is further extended 
when additional point defect scatterers are introduced 
via Ge-alloying. The resulting Sii_a;Gex nanocompos- 
ite material allows both grain boundary and point defect 
phonon scattering mechanisms to be tuned while optimiz- 
ing the efficiency of the resulting fused nanocrystalline 
network. The effect that these added scattering mecha- 
nisms have in the distribution of the phonon density of 
states (PDOS) of the resulting nanocomposites however 
remains an ongoing area of intereslPi'. 

Modifications to the PDOS in nanocrystalline alloys 
due to both finite-size and grain boundary effects have 
been previously explored via a number of probes where 
phon on lifetime broadening effectfl^ as well as both red- 
l ^^ l ^^ l and blue-shift^i^ in phonon frequencies have been 



reported. Additionally, strain fields introduced via the 
crystallite reduction process and frozen at grain bound- 
aries can modify the resulting PDOS through renormal- 
ized force constants due to bond length changes near 
interfaceJi^. Enhancements in the low energy spectral 
weight of the PDOS in a number of systems are often 
ascribed as interface modes or intercrystallite vibrations, 
and lifetime broadening effects are typically observed to 
reduce the sound velocityin the low energy acoustic 
regime of nanocrystalliteaI^"l21. The interplay between 
these and confinement effects associated within the lat- 
tice dynamics of nanocrystals however remains an active 
area of investigation where experimental studies map- 
ping the phonon mode distributions in nanocomposites 
with anomalous thermal transport properties are rather 
scarce. 

In this article, we report inelastic neutron scatter- 
ing measurements mapping the phonon density of states 
across a series of nanocomposite Sii_a;Gea; [x — 0, 0.05, 
0.10, 0.20) samples clustered near the Si-rich end of the 
phase diagram that were previously shown to possess sub- 
stantially reduced thermal conductivitj^. The effect of 
adding point defect Ge scattering sites (in the low concen- 
tration limit) into the hot pressed nanocrystalline matrix 
as well as modifications to the PDOS imparted via the 
nanotexturing process are explored. At higher Ge con- 
centrations, the shift in spectral weight observed upon 
the introduction of Ge into the Sii-^^Ge^ matrix devi- 
ates from expectations of theoretical estimates assuming 
a homogeneous Si-Ge mixing that smoothly renormal- 
izes the lattice force constants. Instead, signatures of 
a two-mode behavior and Ge clustering appear in the 
resulting spectra^-^ as Ge moves far away from the point- 
defect regime. In comparing the PDOS of bulk, poly- 
crystalline Si with hot pressed, nanostructured Si (whose 
thermal conductivity has been reduced by nearly an order 
of magnitude) , the phonon spectra of the two compounds 
appear nearly identical. This suggests a minimal role of 
modified single phonon dynamics above our experimental 
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resolution in changing the heat conduction mechanism of 
this thermoelectric system. 

Our results also allow a comparison between bulk poly- 
crystalline Si and as-milled Si nanopowder at three dif- 
ferent temperatures (300 K, 500 K, and 700 K). While 
the PDOS of bulk Si is nearly temperature indepen- 
dent up to 650 K, an anomalous enhancement in the 
lattice dynamics of the nanopowder Si appears near 35 
meV upon warming. At 300 K, a lifetime broadening of 
phonon modes in the transverse acoustic energy regime 
appears in the nanopowder spectrum, resulting in an en- 
hancement of the low energy PDOS relative to the bulk. 
More generally, the competition between energy scales 
in nanocomposites such as energy raising surface effects 
versus free energy lowering enhanced entropies intrinsic 
to the nanostructures themselves (such as vibrational en- 
tropy, Svi^^ remains unexplored in this material. Our 
measurements allow the calculation of the vibrational en- 
tropy change of Si upon transitioning from the bulk crys- 
talline phase to the highly strained ball milled nanopow- 
der phase and show a sizeable increase in Syib- 



II. THEORETICAL MODEL OF PDOS 

As a guide for interpreting our experimental results, 
the one-phonon density of states (DOS), g(w), is calu- 
clated as: g{uj) — J2j I d{u} — ujj{q))dq where 0Jj{q) 
is the frequency of a phonon mode in branch j with wave 
vector, q, and the integral, taken over the first Brillouin 
zone, is evaluated using the tetrahedron approach. The 
phonon frequencies are obtained by diagonalizing the dy- 
namical matrix. The harmonic interatomic force con- 
stants for Si and Ge are calculated from first principles 
within the framework of density functional perturbation 
theory (DFT) . Here we have used the Quantum Espresso 
progranPil, An 8x8x8 Monkhorst-Pack^^ mesh has 
been used for the electronic states, and a cut-off energy 
of 30 Ry was chosen. The pseudopotentials were gen- 
erated following the method of von Barth and CaP^. 
These pseudopotentials accurately reproduce the mea- 
sured phonon dispersion data for both Si and Ge. The 
exchange-correlation energy is calculated within the lo- 
cal density approximation, using the results of Ceperley 
and Aider as parameterized by Perdew and Zungep^S. 
For the Sii_a;Gea; alloy, we have adopted a virtual crystal 
approximation model, where the interatomic force con- 
stants and atomic masses of pure Si and Ge crystals are 
averaged according to their relative concentrations in the 
alloy. Further details of this approach have been previ- 
ously described elsewher^26H28|_ 



III. EXPERIMENTAL DETAILS 

Three nanocomposite Sii-^^Ge^^ samples with x =0.05, 
0.10, 0.20 were prepared via high energy ball milling 



and hot pressi ng techniques identical to the processes de- 
scribed earlier^. All nanopowder samples were ground 
in a SPEX 8000 high-energy ball mill inside of stainless 
steel jars. These samples were then hotpressed inside a 
custom built press with DC graphite heating elements 
and densified as described earlier^l. We note here that 
each of these hot pressed concentrations had 2% P added 
to the matrix to enable sample compaction and increase 
the carrier concentraion. Each alloy was subsequently 
verified to be single phase via x-ray diffraction measure- 
ments, and EDS measurements confirmed the absense of 
contaminants from the ball milling process to within 1%. 
We note here that hotpressed samples prepared in this 
way are very stable. Previous measurements failed to 
observe significant property and structure changes fol- 
lowing sample annealing at 1100 C in air for over one 
month. 

Our study also looked at three Si samples: (1) a bulk 
Si sample consisting of mortar crushed (5N) Si ingots, 
(2) a nanopowder Si sample created by ball milling the 
same bulk Si ingots in an argon atmosphere, and (3) a 
hotpressed Si sample consisting of Si and 2%P ball milled 
together and then hotpressed in a manner identical to the 
Sii_a;Gea; nanocomposite samples. For the remainder of 
the paper's discussion, we refer to the ball milled and 
hot pressed SiPo.o2 sample as "Hot pressed Si" , the bulk 
mortar crushed polycrystalline Si ingots as "Bulk Si", 
and the loose ball milled Si nanopowder as "Nanopowder 
Si" . All samples were stored under an argon atmosphere 
and eventually loaded into sealed aluminum cans within 
a furnace for neutron measurements. 

X-ray measurements were collected in a Bruker D2- 
Phaser system, and crystallite imaging was performed 
on a JEOL 6340F scanning electron microscope (SEM). 
Neutron measurements were performed on the ARCS 
time-of-flight chopper spectrometer at the Spallation 
Neutron Source at Oak Ridge National LalJ^Sl. A series of 
different incident energies and chopper frequencies were 
utilized to explore the lattice dynamics with variable in- 
strument resolution across different energy ranges. The 
two configurations chosen to map the phonon density of 
states were Ei ~ 120 meV and Ei = 40 meV with in- 
cident Fermi chopper frequencies of 600 Hz and 300 Hz 
respectively. Detector efficiencies were corrected via a 
standard vanadium normalization, and the background 
scattering for a given incident energy and temperature 
was removed by subtracting off from an empty sample 
can under the sample experimental conditions. Samples 
were loaded into cylindrical 1/2 inch inner diameter Al 
cans inside a He atmosphere and sealed prior to measure- 
ment. 

The background subtracted data was summed over all 
Q = - 141-1 ^ -^20 meV) and Q = - SA-^ 
{Ei — 40 meV) for neutron energy loss spectra. Data 
neighboring detector gaps or edges in (Q,E)-space were 
excluded from the summations in order to minimize de- 
tector noise effects within the spectra. Even though Si 
possesses a strong coherent scattering cross section, the 



broad momentum averaging of intensities over multiple 
Brillouin zones generates a reasonable approximation of 
the incoherent phonon density of states, which is verified 
in Section IV of this paper. 

The integrated intensity spectra were then fit via a 
routine similar to that detailed by Kresch et al^. (1) 
removing an initial, constant, multiphonon sum from 
the scattering data, (2) calculating the resulting single- 
phonon scattering profile, (3) determining the corre- 
sponding multiphonon contribution up to 4-phonon scat- 
tering processes by convolving the single-phonon scatter- 
ing profile, (4) scaling this multiphonon contribution up- 
ward to approximate multiple scattering contributions, 
and (5) iteratively repeating this process until conver- 
gence was reached. The multiple scattering contribution 
(which is geometry and packing density dependent) was 
typically of the order unity with the multiphonon contri- 
bution and varies no more than w 7% over the range of 
all samples. For Ei — 120 meV data, a small constant 
background term had to be added to the data due to fast 
neutron/gamma background scattering from the neutron 
source. For all samples, a hard single-phonon cutoff fre- 
quency of 66 meV was enforced, which is supported by 
the raw scattering data. 

In order to compare data of different Ei 's and resolu- 
tion profiles, density of states data from Ei = 120 meV 
and Ei = 40 meV were normalized together at AE = 28 
meV (near a minimum in the PDOS of Si). Specifically, 
in combining the data from each Ei, the total spectral 
weight of the single phonon spectrum in the Ei = 40 
meV data was normalized to the percentage of the den- 
sity of states below 28 meV in the Ei — 120 meV spectra. 
Combing data in this fashion necessitates the approxima- 
tion that the multiphonon background can be smoothly 
scaled; however we stress here that our analysis holds 
for comparing each data set individually as well as com- 
bined. All relative changes reported between each den- 
sity of states spectra are determined within an individual 
Ei configuration, and we cross-checked the normalization 
procedure through comparison of the bulk Si PDOS with 
the DFT model (discussed further in Section IV). 

The resulting generalized phonon density of states 
(GPDOS) is a neutron weighted average over the atomic 
species, d, where the GDPOS cx ^d9di^)ld^ (neglect- 
ing the variation in Debye- Waller factors). Here, gd{E) is 
the partial density of states for atom d, jd is the fraction 
of atomic occupation, (Jd is the total neutron scattering 
cross-section and rud is the mass of the d*'' atomic species. 
For pure Si samples studies, the GPDOS is identical to 
the conventional PDOS. There is an approximately 35% 
decrease in this weighting factor when comparing Ge to 
Si; however the GPDOS derived from our measurements 
of Sii_a;Gea; samples with x < 0.20 remain largely domi- 
nated by the Si scattering. 
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FIG. 1. (a) SEM image of a portion of the Si nanopowder 
sample used in our studies, (b) X-ray diffraction profile from 
the same Si nanopowder sample showing coherent scattering 
domains of approximately 17 nm in diameter. 



IV. GPDOS IN BULK CRYSTALLINE Si AND 
NANOPOWDER Si 

As an initial investigation into the effect of high-energy 
ball milling bulk crystalline Si into a nanopowder, both 
the bulk Si powder and ball milled Si nanopowder were 
explored. SEM images in Fig. 1 (a) show the resulting 
particle size of the Si nanopowder to be « 100 nm; how- 
ever the x-ray scattering profile shown in Fig. 1 (b) shows 
a much smaller correlation length due to a high defect 
density and enhanced strain field imparted via the ball 
milling process. Fits to a Williamson-Hall style plolP^ 
render a coherent diffraction domain size of w 17 nm — 
suggesting a large density of grains within the mean par- 
ticle size resolved via SEM images. Within resolution, no 
volume fraction of amorphous Si was observed in diffrac- 
tion measurements of this or any of the other nanotex- 
tured samples studied in this paper. 

Looking first at the comparative lattice dynamics at 
T — 300 K, the Al-can background subtracted intensity 
maps of both the bulk polycrystalline and nanopowder 
Si samples are plotted in Fig. 2. From the Ei = 120 meV 
data, both samples show a well defined single phonon cut- 
off energy of 66 meV above the transverse optical zone 
boundary and Van Hove singularitj^^HSIi with a broad 
multiphonon peak appearing at ~ 78 meV. From an ini- 
tial inspection of both the Ei = 120 meV and Ei — 40 
meV maps, an overall broadening of spectral features in 
both momentum and energy appears in the nanopowder 
sample relative to the bulk — consistent with the short- 
ened correlation length and lifetime effects of excitations 



4 




02468 10 12 14 2 4 6 8 

Q(A->) 0(A-') 



FIG. 2. (Q, uj) intensity maps showing background sub- 
tracted inelastic neutron scattering spectra with incident neu- 
tron energies Ei = 120 meV for both the (a) bulk Si and (b) 
nanopowder Si samples. Intensity maps with Ei = 40 meV 
are also shown in panels (c) and (d) for the bulk and nanopow- 
der Si samples respectively. 



expected in a highly disordered lattice. In order to quan- 
tify the resulting effect on the GPDOS, the scattering for 
both samples was integrated over momentum and plotted 
in Fig. 3. 

After following the iterative fitting procedure for the 
removal of multiphonon and multiple scattering effects 
detailed in Section II, the resulting singlc-phonon and 
scaled multiphonon components of the total scattering 
are overplotted with the momentum integrated data in 
Fig. 3. The solid line shows the combined scattering 
terms and the resulting fit to the total integrated scat- 
tering. The refinement of the multiphonon scattering 
components is governed by scattering above the single- 
phonon cutoff where a pronounced multiphonon peak in 
the momentum integrated data appears in all data sets. 
This provides a fortuitous guideline for removing multi- 
phonon contamination, and the resulting fits to the data 
model the A_E > 66 meV regime and the multiphonon 
peak near 78 meV well. 

The only exception to this is a high frequency upturn 
and an enhanced constant background scattering term 
observed only in the nanopowder Si sample at 300 K. 
While the enhanced energy-independent scattering term 
can be removed, our multiphonon fits fail to model the 
slight intensity upturn observed for Ai? > 100 meV. Both 
of these spurious high frequency features however disap- 
pear upon heating the Si nanopowder sample, suggesting 
either the influence of an unknown contaminant adsorbed 
to the nanopowder (despite the precaution of maintain- 
ing an Ar environment during preparation and storage) 
or the influence of anomalous grain boundary/strain ef- 
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FIG. 3. Momentum integrated and averaged intensities (black 
squares) for both (a) bulk Si and (b) nanopowder Si sam- 
ples. Closed circles show the reduced single-phonon com- 
ponent of the scattering, and open squares show the mul- 
tiphonon/multiple scattering component of the total inten- 
sities. Solid lines show the combined single phonon, scaled 
multiphonon, and energy independent background intensities 
fit to the total scattering data. 



fects in the phonon spectrum. 

The resulting PDOS obtained from the single-phonon 
profile from the bulk Si sample is plotted in Fig. 4 (a). 
As described in Sect. II, data from the Ei = 120 meV 
dataset transitions to the Ei — 40 meV dataset below 
AE' = 28 meV (denoted by a dashed line) in order to uti- 
lize the superior resolution afforded by a lower incident 
energy over the acoustic phonon regime. As a reference 
we have overplotted the results of our DFT calculations 
(convolved with the instrument resolution) for the PDOS 
of Si with the experimentally determined PDOS. Aside 
from a small frequency shift due to the systematic under- 
estimation of zone boundary energies in the DFT calcu- 
lation, the experimental PDOS matches remarkably well 
with the theoretical expectation. This provides an added 
degree of confidence in the experimental data reduction 
and normalization procedure as well as the momentum 
summation's incoherent approximation to the PDOS. As 
a further check, the low energy PDOS is plotted versus 
E^ in the inset of Fig. 4(a). In the long wavelength 
limit, the PDOSoc E'^ as expected; however this approx- 
imation breaks down quickly above IS.E = 9 meV in the 
experimental data and above /S.E « 5 meV in the DFT 
simulation. 

Now applying an identical reduction procedure to both 
bulk and nanopowder Si samples at higher temperatures, 
the temperature dependence of the PDOS for each sam- 
ple at T = 300 K, 500 K and 700 K are plotted in 
Figs. 4 (b) and (c). As expected, the resulting PDOS 
for bulk Si remains largely unchanged upon heating, al- 
though a slight anharmonicity and redshift of the opti- 
cal zone boundary zone boundary Van Hove feature ap- 
pears at 650 K. In contrast, the nanopowder Si sample 
shows more dramatic changes in the PDOS upon heat- 
ing. At increased temperatures, an anharmonic redshift 
again appears in the optical zone boundary feature of the 
nanopowder; however there also appears a redistribution 
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FIG. 4. (a) Black circles show the reduced PDOS for bulk Si 
at 300 K. The red solid line shows the DFT calculated PDOS 
for Si as described in the text. Inset shows the low energy 
PDOS as a function of where symbols show the PDOS of 
Si and the solid black line shows a linear fit to the low energy 
PDOS. (b) Reduced PDOS for bulk Si are compared at 300 K, 
500K, and 650 K. (c) The reduced PDOS for nanopowder Si 
are compared at 300K, 500K, and 700 K. Dashed line denotes 
the energy crossover between Ei = 120 meV and Ei — 40 
meV data, (d) Q=(l,l,l) Bragg peak widths plotted as a 
function of sample temperature for the nanopowder and bulk 
Si samples (errors within symbols). 
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FIG. 5. Expanded comparison of the 300 K PDOS of bulk 
(open squares) and nanopowder (solid circles) Si below 30 
meV. The solid line shows the measured PDOS of bulk Si 
convolved with a damped harmonic oscillator lineshape as 
described in the text. 



of spectral weight resulting in an enhanced PDOS for 
30 < < 50 meV. The PDOS of the nanopowder sample 
in the acoustic regime, while broadened at 300 K, sharp- 
ens slightly upon heating likely due to grain growth and 
annealed strain effects. 

We note here that the sample's temperature history 
went from 300 K to 700 K and finally back down to 500 
K. This suggests that the increased sharpening of the 
acoustic mode distribution from the 700 K and 500 K 
data is likely due to the continued annealing after the 
700 K data collection (« 4 hrs). As a rough check, the 
peak widths of the Si Q=( 1,1,1) Bragg reflections were 
analyzed for both the nanopowder and bulk Si samples 
as a function of temperature and plotted in Fig. 4 (d). A 
slight narrowing of the nanopowder Bragg widths was ob- 
served upon heating to 700 K which supports the notion 
of the acoustic mode lifetime sharpening as arising from 
an annealing/grain growth effect. The peak widths how- 
ever remain notably broader than those of bulk Si (the 
spectrometer's resolution) after heating, demonstrating 
that the sample remains in a largely nanotextured phase. 

The spectral weight enhancement above 30 meV in the 
nanopowder however is robust and scales monotonically 
with temperature, suggesting an origin intrinsic to the 
nanostructured sample. This contrasts with the naive 
expectation of an annealing effect such as that observed 
at lower energies which should smoothly transition the 
dynamics from the nanopowder to those of the bulk. At 
present, the origin of this enhancement is unclear, al- 
though grain boundary modes or defects may play a role. 
Specifically, these modes can in principle be unpinned 
and thermally activated upon heating and would renor- 
malize the spectral weight of the PDOS distribution if 
the volume fraction of defects was large enough. 
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FIG. 6. Ei — 120 meV data showing the resulting single 
phonon, scaled multiphonon, and total scattering fits over- 
plotted with the momentum averaged intensity data for (a) 
Sio.95Geo.05, (b) Sio.9oGeo.io, and (c) Sio.80Geo.20- As a refer- 
ence, overplotted in each panel are the corresponding single 
phonon, scaled multiphonon, and momentum averaged total 
scattering intensities for the Ge-free hot pressed Si sample. 



One striking feature in comparing the PDOS of the 
nanopowder and bulk Si samples is the preferential 
broadening of the spectral weight distribution in the 
transverse acoustic regime. There is little difference in 
the spectral shape near the Van Hove singularities in the 
optical regime of the PDOS, but a clear phonon broad- 
ening is resolvable below 30 meV. The broadening at low 



energies is highlighted in a comparison of the bulk and 
nanopowder samples below 30 meV plotted in Fig. 5. In 
order to estimate the magnitude of broadening, a conven- 
tional method of modeling the spectrum is to convolve 
the ide al spe ctrum with a damped harmonic oscillator 
functiorpSED, The most practical approach is then to 
first consider the measured bulk Si phonon PDOS as a 
model starting point. 

We have already shown that the calculated PDOS con- 
volved with the instrument resolution matches the mea- 
sured spectrum for bulk Si well. Therefore, we can as- 
sume the intrinsic line widths of the bulk Si spectrum 
are well below the instrument resolution, and we can 
convolve the function T{E) cx ^E^^E-E/kr+i^/W' ^ ^'^ ^ ^^ ^ 
with the experimentally obtained Si PDOS to determine 
the needed /? to approximate the nanopowder response. 
The best fit obtained for the nanopowder yields a result- 
ing /3 = 35. The quality factor, /3 is proportional to the 
number of wavelengths a phonon of a given frequency is 
able to propagate within a crystallite domain, and the 
relatively high value for our Si nanopowder can be at- 
tributed to the l arger Si crystallite sizes relative to those 
in earlier studieJ^^EZl. 



V. VIBRATIONAL ENTROPY CHANGE IN 
NANOPOWDER Si 

The vibrational entropy change is definecP^ as AS — 
-3kB {PDOSnano- PDOSbuik)ln{E)dE. In looking 
at the reduced PDOS for Si upon transitioning from the 
bulk to the nanopowder phase, a /S.S^qqk = 0.84 ± 0.05 
ks/atom is calculated from the data at 300 K. This un- 
usually large entropy difference is likely an artifact from 
a slight over subtraction of the anomalous 300 K back- 
ground scattering at low energies. If we instead look at 
only the phonon spectra derived from the Ei — 120 meV 
data above AE = 8 meV, then AS300K = 0.21 ± 0.03 
kB/atom — in line with value s observed in a number of 
metallic nanocomposite^^^l^. Upon warming to 700 K, 
the anomalous background term in the nanocrystalline 
sample's scattering profile vanishes, and the AS^qk re- 
mains nearly unchanged at 0.20 ± 0.03 kB/atom. After 
then cooling down and measuring at 500 K, the resulting 
AScjQQx retains a similar value of 0.23 ± 0.03 kB/atom. 

The differences in the lattice heat capacities be- 
tween the bulk and nanopowder Si can also be cal- 
culated at each temperature through the relation: 
a ^ 3i? PDOS{u:){f^f (//iXip du [J mole-' 
K^']. As an initial check, the heat capacity of bulk Si at 
300 K can be calculated using the measured PDOS giv- 
ing Cy^buik = 20.3±0.4. This agrees reasonably well with 
known value of C„ = 19.9 determined from calorimetry 
measurement^^. At 650 K, the heat capacity for bulk 
Si increases to Cv_buik = 24.1 ± 0.5 (near the expected 
value of 24.^^ while for the nanocrystalline sample the 
heat capacity shows a slightly smaller value of C'y_nano = 
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FIG. 7. Reduced 300 K GPDOS for the hot pressed (a) 
Sio.95Geo.05, (b) Sio.9oGeo.io, and (c) Sio.80Geo.20 samples. 
Solid lines are the corresponding DFT model calculations 
for the PDOS assuming the virtual crystal approximation. 
Dashed line again denotes the crossover from Ei = 120 meV 
data to Ei = 40 meV data. 



23.9 ± 0.4. After cooling to 500 K, C^^buik = 23.5 ± 0.4 
(23.6 expected^^) and Cy^nano = 23.1 ± 0.3. Naively, 
one would expect the heat capacity of the nanopowder 
Si to be higher than that of the bulk Si based on the non- 
neglible difference in the vibrational entropy between the 
two; however the slight decrease in Cy arises due to the 
preferential weighting of higher frequency differences in 
the PDOS by the ln{E) term in the AS calculation rel- 



ative to the Einstein distribution weighting term in the 
Cv integral. 



VI. GPDOS IN NANOSTRUCTURED Sii^Ge^ 

We now turn to the evolution of the GPDOS in hot 
pressed, nanostructured thermoelectric Sii_a;Ge^ alloys 
with a; = 0, 0.05, 0.10, 0.20. The same data collection 
and reduction procedure previously described for bulk 
and nanopowder Si samples was employed on these hot 
pressed thermoelectric samples. The resulting momen- 
tum integrated total scattering plus the single- and mul- 
tiphonon components for each sample at T = 300 K are 
plotted in Fig. 6. Here each Ge-doped concentration is 
overplotted with the scattering components from the hot 
pressed Si sample as a reference for changes in the raw 
scattering profile. Solid lines are again the composite 
fit to the total scattering, and the scaled multiphonon 
background fits the scattering beyond the single phonon 
cutoff well in all cases. The multiphonon contribution 
changes only slightly upon adding Ge to the system, and 
the largest changes occur at low energies as the single 
phonon spectral weight is redshifted downward. 

The resulting, reduced GPDOS spectra for each Ge- 
doped sample are plotted in Fig. 7. Overplotted with 
each GPDOS spectrum are the results of DFT calcula- 
tions of the PDOS for the corresponding Sii-^jGe^, alloy 
assuming a uniform Si and Ge distribution. As Ge is 
progressively alloyed into the Sii_a;Ge^ matrix, the spec- 
tral weight near the longitudinal acoustic zone bound- 
ary is enhanced between 30 and 40 meV. At lower ener- 
gies, there appears an overall broadening of the acoustic 
phonon mode distribution relative to the parent Si mate- 
rial for small Ge concentrations {x = 0.05 and x = 0.10) 
and an overall shift to lower frequencies for the x = 0.20 
sample. The resulting acoustic mode distribution for the 
X = 0.20 composition suggests a lifetime broadened form 
of the DFT theoretical expectation; however a simple 
convolution of the DFT modeled PDOS with a damped 
harmonic oscillator line shape fails to model the data. 

For X < 0.10, the response of the GPDOS to the addi- 
tion of Ge mimics a simple lifetime broadening term in 
the PDOS of Si with no appreciable frequency shift with 
the exception of a low energy shoulder that builds with 
increased Ge concentration. We note here that the nor- 
malized presentation of both Ei = 40 and Ei = 120 meV 
data fails at the 28 meV transition region for the x=0.20 
sample due to appreciable spectral weight present in- 
duced at the crossover point; however the above analysis 
holds since relative changes are only ascertained within 
a single Ei setting for a given sample. 

Now turning to the higher frequency optical phonon 
modes, the addition of Ge smoothly suppresses the spec- 
tral weight of the 60 meV zone boundary modes corre- 
sponding to Si-Si vibrations. The characteristic energy of 
the optical Van Hove singularity is left unchanged how- 
ever. This implies substantial clustering effects in Ge- 
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FIG. 8. (a) Comparison of the 300 K PDOS of the three pure 
Si samples in bulk (square symbols), hot pressed (diamond 
symbols), and nanopowder (circle symbols) forms, (b) Over- 
plot showing the comparison of the 300 K GPDOS for the 
series of hot pressed thermoelectric Sii-a;Gei, samples with 
x=0, 0.05, 0.10, and 0.20. Dashed line again denotes the 
crossover from Ei — 120 meV data to Ei — 40 meV data. 



doped samples and a resulting two-mode behavioi'^SlEIl, 
illustrated most clearly by the appearance of a second 
optical phonon peak in the GPDOS at 50 meV m the 
X = 0.20 sample. In all samples however, there remain 
substantial regions of Si clusters that retain their char- 
acteristic Si-Si 60 meV peak. 



lower energy cutoff of AE = 4 meV or possess relative 
spectral weights below our experiments' sensitivity. 

In comparing the nanopowder Si sample to the bulk 
and hot pressed Si samples, there appears a substan- 
tial lifetime broadening effect resolvable in the acoustic 
regime which results in an enhanced tail in the PDOS 
extending down to the lowest energies measured. The 
effects of this lifetime broadening are removed once the 
nanopowder is hot pressed into a functional alloy; how- 
ever the nanopowder PDOS spectrum may serve as a use- 
ful tool for amplifying otherwise subtle effects imparted 
via the high-energy ball milling process. 

Now turning to the 300 K GPDOS data for the hot 
pressed Sii-ajGe^; alloys, a clear two-mode behavior ap- 
pears as Ge is progressively doped into the materiaP^l. 
Phonon modes at the Si-Si transverse optical zone bound- 
ary retain their characteristic frequency while their rel- 
ative spectral weight is damped with increasing Ge con- 
tent. This suppression in the spectral weight around sa 60 
meV is compensated by a broad redistribution/smearing 
of weight to frequencies as low as 30 meV. At energies 
below 25 meV, the spectral weight distribution largely 
attributable to the transverse acoustic zone boundary 
modes in pure Si broadens with small Ge concentrations 
before undergoing a global shift to lower frequencies at 
X = 0.20. This implies that the clustering of Ge within 
the matrix at low concentrations primarily serves as a 
point defect scattering mechanism for reducing acous- 
tic phonon lifetimes. At large enough densities of Ge- 
impurities however, the lattice undergoes a global re- 
sponse and dispersion relations are renormalized. 

The experimental picture of point defect scattering 
versus grain boundary scattering effects in Sii_a;Gea; 
therefore suggests that the Ge point defect scatterers, 
even at concentrations as low as 5%, introduce a sub- 
stantial change to the global GPDOS relative to simply 
hot pressing a frozen network of grain boundaries in pure 
Si. At low enough Ge concentrations, the overall effect on 
the GPDOS at acoustic phonon frequencies empirically 
mimics that observed in pure nanopowder Si; however the 
potential thermal activation of grain boundary modes or 
defects milled into the Si nanocrystalline matrix at high 
temperatures and frequencies remains an issue for future 
experiments. 



VII. DISCUSSION 



VIII. CONCLUSIONS 



A summarized comparison of the GPDOS determined 
for each sample measured at 300 K is show in Fig. 8. 
In comparing the PDOS of the bulk, hot pressed, and 
nanopowder Si samples in Fig. 8 (a) there is no dis- 
cernible shift in the spectral weight distribution of the 
PDOS between the bulk and hot pressed Si samples. This 
stands in stark contrast to dramatic reduction in the ther- 
mal conductivity observed in the hot pressed Si sample^, 
implying boundary modes and other scattering effects 
must take place at energies well below our experimental 



To conclude, we have investigated the evolution of the 
GPDOS in a series of nanostructured Sii-xGe^ thermo- 
electric alloys as well as changes to the PDOS within 
nanopowder Si created via high-energy ball milling. Ge- 
substitution into Sii_a;Gea; nanostructured alloys results 
in a two-mode type of behavior where clustering effects 
are evident. For light Ge-concentrations with x < 0.10, 
the acoustic phonon mode distributions broaden consis- 
tent with a substantial reduction in the acoustic phonon 
lifetime; however hot pressed, nanostructured Si exhibits 
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a nearly identical PDOS to that of bulk Si despite an or- 
der of magnitude reduction in thermal conductivity. This 
implies the dominant role of long wavelength and rare- 
region scattering mechanisms beyond the resolution of 
our experiments. Our experiments do however resolve a 
dramatic lifetime broadening of acoustic phonon modes 
in pure nanopowder Si relative to the bulk along with an 
increase in the vibrational entropy of AS* w 0.20 ± 0.03 
kB/atom. A temperature dependence was reported in 
the PDOS of nanopowder Si distinct from the bulk sys- 
tem that likely implies the thermal activation of addi- 



tional modes within the nanocrystallite lattice. 
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